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A quantum mechanical model for the description of the a-decay of heavy nuclei with 
accompanying photons emission is presented. The model is based on a quantum mechan- 
ical one-particle model of a-decay through a decay barrier. The bremsstrahlung spectrum 
calculation employs multipole expansion of the vector potential of the Coulomb field of the 
daughter nucleus and takes into account the dependence on the angle between the directions 
of the a-particle propagation and the photon emission. Spectra of 210 Po are obtained for 
the angles 25° and 90°, and the best agreement with the experimental data in the 90° case 
^\ ' in a comparison with other existing models is achieved. From the angular analysis, the 

' model gives monotonic behavior of the bremsstrahlung spectrum for any value of the angle. 

We find that sub-barrier photon emission (i. e. bremsstrahlung photon emission during a- 
particle tunneling through decay barrier) exists but gives a small contribution to the total 
^ ' bremsstrahlung spectrum. 

o : 
o\ ■ 

O ■ §1. Introduction 

The existence of sub-barrier bremsstrahlung phenomena in a-decay (i.e., brems 
-strahlung photon emission during tunneling of an a-particle through a decay barrier) 
has not yet been confirmed through the analysis of experimental data. In Ref. 1) 
' we proposed a quantum mechanical approach to resolve this problem. With this 

approach, one can calculate the bremsstrahlung spectra both with and without sub- 
barrier photon emission for the nuclei 210 Po, 214 Po, 226 Ra and 244 Cm, for which 
experimental data 2 ^ 5 ^ exist, and with a comparative analysis determine whether 
J> ■ sub-barrier bremsstrahlung exists. 

This approach can be realized only with such a model that yields sufficiently 
good agreement with experimental data. Note, that the present theoretical descrip- 
tion of experimental data 2 ) for 210 Po is unsatisfactory, and there is a difficulty in 
explaining the existence of the "hole" in the bremsstrahlung spectrum. Note that 
these experimental spectrum carried out for a value of the angle between the direc- 
tion of the photon emission and the direction of the a-particle propagation equal 
to 25°. New experimental data 6 ) of the bremsstrahlung spectrum for 210 Po have 
been obtained for the different value of the angle equal to 90°, exhibit monotonic 
behaviour and have no "hole" with such amplitude. We suppose that obtaining a 
description of the bremsstrahlung spectrum for 210 Po with a sufficient degree of ac- 
curacy and clarification of the spectrum behaviour are important problems in the 
investigation of sub-barrier bremsstrahlung phenomena. 

To make possible a comparative analysis of these two groups of experimental data 
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for °Po obtained for different angular values, we present here an angular quantum 
mechanical model of the a-decay of heavy nuclei with accompanying (spontaneous) 
emission of photons. The foundation of this work was laid in Ref. 7), where we 
calculated the bremsstrahlung spectrum for the nucleus 214 Po, taking into account 
El. 

This model is based on a quantum mechanical one-particle model of a-decay 
through a decay barrier. 8 ) A description of the electro-magnetic field of the daughter 
nucleus and the emission of photons in the a-decay is constructed on the basis of 
quantum electrodynamics. 9 ) In calculating the bremsstrahlung spectrum, we use 
an expansion of the vector potential of the electro-magnetic field of the daughter 
nucleus in electric and magnetic multipoles, 10 ) in contrast to Papenbrock and Bertsch 
quantum mechanical method, n ) and we do not use the additional transformation 
(2- 14) used in Tkalya model. 12 -* We first present the mathematical tools of this 
model. 

§2. A non-stationary quantum mechanical model of a-decay 

with bremsstrahlung 

For investigation of the photon emission process, we consider the state of the 
system (the a-particle and the daughter nucleus) before the emission as the initial 
(i-state) and the state of the system after the emission as the final (/-state). We 
calculate the bremsstrahlung spectrum using the transition matrix element 13 ) 

t 

a fi = (k f ,n k + l\V(r, t)\ki,n k ), V(r,t) = -i Je iAot 'v{r,t')e- iAot 'dt', (2-1) 

o 

where tpi(r) = \k^ and ipf( r ) = \kf) are the wave functions of the system in the 
initial i-state and the final /-state from the unperturbed operator Hq which describes 
scattering of the a-particle upon the daughter nucleus, uses the barrier and does not 
take into account the photon emission; w = E; n k is the number of photons of one 
sort with impulse k in the initial z-state. We choose units for which h = 1, and 
c = 1. In the Coulomb calibration, the interaction operator has the form 10 )' 13 ) 

V(r,t)=-Z eS -Ap, (2-2) 
m 

where Z e R is the effective charge, m is the reduced mass of the system, and A is the 
vector potential of the electromagnetic field of the daughter nucleus. 
We choose the vector potential A(r,t) to take the form 9 ) 

A(r,t) = ^(c k , a A k!a + c+ a At ta ), A k<a = JfeW^M, (2-3) 

k,a V / \ W 

where e^ a ) is the unit polarization vector of the photon, k is the wave vector of the 
photon, and w = k = |fc|. The vector e( Q ) is perpendicular to k in the Coulomb 
calibration. We have two independent polarizations, e^ 1 ) and e^ 2 ) for any photon 
with impulse k. 
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Let us consider a wave packet of the form 

+00 

*i,/(r, t) = J g(k - hj)il>i,f(r, tye'^dk, (2-4) 


where in the expressions ^,/(r, t) and k^j we use the index % or / only in dependence 
on consideration of the initial %- or final /-state of the system. We take the weight 
amplitude g(k — kij) in the form of a Gaussian: 14 ) 

g(k-k iJ ) = Ce W (- { \^kf ) (2 " 5) 

where C is a normalization factor and can be obtained from the normalization con- 
dition J \g(k - kij)\ 2 dk = 1. 

The wave packets (2-4) can be used for analyzing in time the tunneling of the 
a-particle through the decay barrier and allow to calculate the tunneling time of this 
particle in dependence on the energy level of such process. Preliminary calculations 
of these times of the a-decay with taking into account a possibility to emit photons 
from the barrier region are presented in Ref. 7). Development of the model of the 
bremsstrahlung in the a-decay on the basis of the wave functions in form (2-4) allows 
to calculate both the bremsstrahlung spectra and tunneling times. Therefore, we use 
the wave packets (2-4) for ipi(r) and ipf(r), with the time factor exp (—iwt). 15 ^ 

Substituting this wave packet into (2T) and taking account of the property 

e - iAot ipi(r) = e- iWit ipi(r), ip* f (r)e iAot = ^* f {r)e iw f\ (2-6) 

we obtain the matrix element in the form 

+00 +00 

dfi = J dki j dk2g(h - ki)g(k 2 - kf)*(k2,n k + l\V(t)\ki,n k ), (2-7) 


where 

t 

V{t) = -i j e iW2t 'Ve- iwit 'dt'. (2-8) 


We study the spontaneous emission of one photon with an impulse k. Therefore, 
rif, = 0, and the sum ^ can be omitted in the calculation of the matrix element (2-1). 

k 

For calculation of the bremsstrahlung spectra we use the stationary approxima- 
tion: t — > +00. Taking into account (2-2), (2-3), (2-7) and (2-8), and also using the 
definition of the 5-function, we calculate the transition matrix element 

(fc2,l|V"|fci,0) = F 21 27r5(w 2 + w - Wl ), (2-9) 

where 

e l2ir 
F 2 i = Z eS — \ — p(k 1 ,k 2 ), 
m V w 
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p(h,k 2 )= V e^>(h 



a=l,2 



—ikr 



d_ 
dr 



For monochromatic particles, we obtain 9 )' 15 ) 

I 1 2 

afi = 2irFf i | C I S(w + Wf — Wi). 



(2-10) 



(2-11) 



We define the transition probability from the initial z-state to the final /-state 
in a unit time with a single photon emission of impulse k and polarization e*") as 
follows: 9 ^ 



mk f w%\F fi \ _ Z 2 s e 2 k f w fl 



dW _ 

~d!T v ~ (27) 3 (2vr) 2 
Wfi = Ei- E f . 



(27r) 4 m 



\p(ki,k f )\ 2 , 



(2-12) 



The most difficult step in calculating the bremsstrahlung spectrum is calculating 
p(ki, kf). We have 

+oo 

p(h,k f )= e(a> f dr f df2r 2 r f (r)e- ikr ^i(r). (2-13) 

a=l,2 o r 

In calculating p(ki,kf), various approaches can be used, and the various quan- 
tum mechanical methods for these approaches differ. Papenbrock and Bertsch used a 
dipole approximation in obtaining a description of the vector potential of the electro- 
magnetic field. 11 ) Tkalya used a multipole approach in calculating the bremsstrahlung 
spectrum, 12 ) taking into account El and E2. However, in calculating the integral 
(2-13) both in Tkalya model and in Papenbrock and Bertsch model the transforma- 
tion 



■(f\p\i) = (f\[H, P ]\i)— = ih(f 

Wfi 



I Wfi 



(2-14) 



is used. This transformation can be used only in a case that the approximation 

exp(ifcr) = 1 (2-15) 

is applied (see Appendix A). Therefore, using the transformation (2-14), the accuracy 
of the calculation of the integral decreases, and the calculations of the spectra with 
the models 11 )' 12 ) are carried out not in the dipole approximation or the multipolar 
approximation (taking into account El and E2) but in the approximation (2-15). 
For this reason, we do not use (2-14). 

In computing the integral (2-13), we apply the multipolar expansion of the vector 
potential of the electromagnetic field without using the transformation (2-14). We 
write the polarization vectors e a in terms of circular polarization vectors £ with 
opposite directions of rotation: 9 ) 



V2 



(2-16) 



Does Sub-Barrier Bremsstrahlung in a-Decay of 210 Po Exist? 



5 



We obtain 



+00 

p(ki,k f )= Y, I dr I 'df2r 2 r f (r)e- ikr ^-Mr), (2-17) 

where 

/i_i = -^(l-i), ^ 1 = _L(_i_i). (248) 
For the initial i- and final /-states, we can write 

A,f(r)=^ij(r)Y lm (K' f ), (2-19) 

where Yi m {n^) are the normalized spherical functions. Following Ref. 10), we have 

d , . . dtbAr) m . ,-. 
^(r) = -^ToloK), 

1 

T i,oK) = 2 ( 110 l " wO)y l5 _ M K)e M , (2-20) 

M=-l 

where are the vector spherical harmonics, and (j a lj| — W*0) are the Clebsch- 

Gordan coefficients. 

We expand the vector potential A(r, t) in multipoles as 10 ) 

i^ kr = v^F/x J] (v / 2TTTi / [A, m (r, M) + ^A, ro (r, £)A , (2-21) 
z=i ^ ' 

with 

A im (r,M) = ji {kr)Tu ^{n r ), 

A lm (r,E) = yj ~^^ii- 1 {kr)T n -i^{n r ) ~ \j ^f—ji+i(kr)T ll+1 ^(n r ), 

(2-22) 

where ji(kr) is the spherical Bessel function of order I, and A; m (r, M) and A; m (r, 22) 
are the magnetic and electrical multipoles, respectively. 
From the calculations, we obtain p(h, kf) as 



p(fe J ,A; / ) = v^^(v^TTHyb Mi (fe J ,fe / )-^(fe i ,A ;/ )] 
l=i ^ 

p M \k u k f ) = hJ l {l), 

p El (k u k f ) = -JB±I 2Jl (l - 1) + ^JZ hJl{l + i), (2 . 23 ) 

where 

+00 

Ji(") = / r 2 ^,(r)^li n (A;r)dr, (2-24) 
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n=-i 





(2-25) 



n=-i 

and n G N (here N is the natural number set), I is the quantum number for the 
wave function of the final /-state. For the transitions Ml and El we obtain 



where Ad and are the mass number and charge of the daughter nucleus, and 6 
is the angle between the direction of the a-particle propagation and the direction of 
the photon emission. 

§3. Angular analysis of the bremsstrahlung spectrum in I^Po a-decay 

We calculated the bremsstrahlung spectra for the a-decay of g4°Po at angular 
values 25° and 90° using the model described above. Our results for 90° are presented 
in Fig. 1. In the calculations, we used the parameter values of the Coulomb barrier 
from Ref. 11) and the wave function of the initial z-state as the decay state and the 
final /-state as the scattered state. 

In Fig. 1, it can be seen that the line 3 obtained using our model for 90° lies 
nearer the experimental data 6 -* than both the line 5 representing the bremsstrahlung 
spectra obtained using the models of Papenbrock and Bertsch 11 ) and Tkalya 12 ) (tak- 
ing into account El), based on quantum electrodynamics, and line 1 representing 
the bremsstrahlung spectrum obtained using Tkalya model 12 ) based on classical 
electrodynamics. Such an accurate description of the experimental data 6 ) of the 
bremsstrahlung spectrum for the nucleus 210 Po has been obtained using a multipole 
approach without the transformation (2-14). 

The total bremsstrahlung spectrum for 210 Po calculated for 90° using our model 
is represented by line 3 in Fig. 1 and is seen to be monotonic. The spectrum calcu- 
lated for 25° lies close to it. With angular analyzing, our model is found to predict 
a monotonic dependence of spectra on photon energy for any angle value (in the 
region 5° - 90°) and does not yield a slope as large as that of the experimental 
data. 2 ) Taking into account Ml, E2 and M2 multipoles the monotonic behaviour 
of the bremsstrahlung spectra for any angular value exists also where there is no 
visible "hole" with such amplitude, as in the experimental data. 2 ) As we have seen, 
the quantum mechanical models of Papenbrock, Bertsch 11 ) and Tkalya 12 ) also fail 
to explain the existence of the "hole" in the spectrum. 




(2-26) 



Z e & = 



2A d - AZ d 
A d + 4 



(2-27) 
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Fig. 1. Bremsstrahlung in the a-decay of g4°Po: 1 — the result of Tkalya model 12 - 1 based on classical 
electrodynamics; 2 — experimental data; 6 - 1 3 — the result of our model for 90° with but El 
and sub-barrier photon emission; 4 — the result of our model for 90° with El but without 
sub-barrier photon emission; 5 — the result of Tkalya quantum mechanical model taking into 
account El 12 - 1 and that of Papenbrock and Bertsch model; 11 - 1 6 — experimental data. 2 ' 



Note that taking account of the mixed region 16 ) is useful for detailed analy- 
sis of photon emission from different regions and the contributions to the total 
bremsstrahlung spectrum. However, we find that the total bremsstrahlung spectrum, 
calculated on the basis of the quantum mechanical models 11 )' 12 ) and our model, is 
not depended on the choice of the mixed region and its change. 

In addition, we calculated the bremsstrahlung spectrum for 210 Po without taking 
into account the sub-barrier emission of the photons. The resulting spectrum is 
represented by line 4 in Fig. 1 for 90°. Here, it can be seen that line 4 lies somewhat 
farther from the experimental data 6 ) (obtained for 90° also) than line 3 for the 
total bremsstrahlung spectrum. From this we can conclude that taking account of 
the photon emission from the barrier region in the a-decay of 210 Po increases the 
accuracy of experimental spectra description. But sub-barrier bremsstrahlung gives 
a small contribution to the total bremsstrahlung spectrum. 

§4. Conclusions and perspectives 



We presented a quantum mechanical approach to the calculation of the brems- 
strahlung spectrum in the a-decay of heavy nuclei, where the angle between the di- 
rections of the a-particle propagation and the photon emission is taken into account. 
On the basis of the spectrum calculation carried out with this approach, we conclude 
that the model gives a monotonic bremsstrahlung spectrum for the different angles. 
(For readers interesting in the comparative analysis of the theoretical or experimen- 
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tal bremsstrahlung spectra obtained for different angles, we can also propose a new 
method presented in 17), where the dependence of the bremsstrahlung spectra on 
the values of such angle is obtained in a simple analytical form and is more obvious 
and simple.) 

According with analysis by our model, taking into account multipoles of larger 
orders does not cause the spectra to become non-monotonic and does not explain 
the existence of the "hole" in the experimental data 2). The "hole" in the total 
bremsstrahlung spectrum reported in Ref. 2) cannot be explained by the choice of 
i?2 ; its adjustment, the arrangement of the mixed region (which was introduced in 
Ref. 16)) nor use of the isotropic quantum mechanical model, where the potential 
starting from Ri is pure Coulomb. (Here, Ri and R2 are the internal and external 
radii of the Coulomb barrier. 7 )' 12 )' 16 )) We agree with Tkalya that it would be useful 
to obtain more detailed experimental data for the bremsstrahlung spectrum in the 
a-decay of 210 Po for different angle values. 

Prom a comparative analysis, it is found that taking into account the sub-barrier 
bremsstrahlung phenomena in a-decay of the nucleus 210 Po increases the accuracy 
of experimental spectra description. But sub-barrier bremsstrahlung gives a small 
contribution to the total bremsstrahlung spectrum (see Fig. 1). From this result, we 
can explain the fact that the models constructed on the basis of classical electrody- 
namics and those employing a semiclassical approximation describe the experimental 
data for the bremsstrahlung spectra of other nuclei sufficiently well. 

We would also like to note further matters regarding the investigation of sub- 
barrier bremsstrahlung phenomena. It would be useful to obtain more accurate 
experimental data for the nucleus 210 Po for photons energies Wfi > 400 MeV and for 
the different angle values. The accuracy of the description of the a-decay with photon 
emission can be improved by accounting for the individual properties of the nucleus. 
According to our model, small oscillations exist in the bremsstrahlung spectrum for 
210 Po and the period can be estimated. The existence of oscillation can be explained 
through use of the multipole approach without recourse to the transformation (2-14) 
in the calculation of the bremsstrahlung spectrum. 
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Appendix A 

Calculation of Multipoles 

In calculating the bremsstrahlung spectrum, we use a matrix element of the form 

(f\pe ikr \i). (A-1) 
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In the calculation of this quantity, Tkalya 12 ) and Papenbrock and Bertsch 11 ) used 
the following transformation: 



-(f\p\i) = (f\[H,p]\i)—=ih{f 
w fi 



'•> — ■ (A-2) 



However, this transformation can be used only in case that approximation 

e ikr = 1 (A-3) 
is made. Indeed, taking into account the properties 

Hf(r)^(r) = (J?- A + U(r)j f{r)^{r) + f(r) (J?- A + U(r)^j Y>(r) = /(r)^(r), 

(A-4) 

we have 

(f\Hpf(r)\i) = {f\E f pf(r)\i) = E f (f\pf(r)\i), 
{f\pHf(r)\i) + (f\pf{r)H\i) = (f\pf(r)E t \i) = E t (f\pf(r)\i). (A-5) 

Only in the case f(r) = const does the inequality (A-4) become an equality, in which 
case the transformation (A-2) can be used to simplify the calculation of (A-l). In a 
multipole expansion, we have 

+°° m=l /L.\ /LA 

f(r) - e ikr = ^Y.YZ il ^ kr ) Y im T ^ U ' (A " 6) 

1=0 m=-l ^ ' ^ ' 

where ji(kr) is the spherically symmetric Bessel function of order /. 
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